Adiabatic Effective Action for Vortices in Neutral and Charged 

Superfluids 



M. Hatsuda 1 , M. Sato 2 , S. Yahikozawa 2 and T. Hatsuda 3 
1 National Laboratory for High Energy Physics, Oho 1-1, Tsukuba, Ibaraki 305, Japan 
2 Physics Department, Kyoto University, Kyoto 606-01, Japan 
3 Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 

Abstract 



Adiabatic effective action for vortices in neutral and charged superfluids at 
zero temperature are calculated using the topological Landau-Ginzburg the- 
ory recently proposed by Hatsuda, Yahikozawa, Ao and Thouless, and vortex 
dynamics are examined. The Berry phase term arising in the effective action 
naturally yields the Magnus force in both neutral and charged superfluids. 
It is shown that in neutral superfluid there is only one degree of freedom, 
namely the center of vorticities, and the vortex energy is proportinal to the 
sum of all vorticities so that it is finite only for the vanishing total vorticity 
of the system. On the other hand the effective mass and the vortex energy 
for a vortex in charged superfluids are defined individually as expected. The 
effects of the vortex core on these quantities are also estimated. The possible 
depinning scenario which is governed by the Magnus force and the inertial 
mass is also discussed. 
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I. INTRODUCTION 



Vortex dynamics play an essential role in various properties of superfluids and supercon- 
ductors Recently, numbers of experiments on the Hall effect [§J and the anomalous 
Hall effect || in superconductors have provided renewed interest in the theoretical study 
of vortex dynamics p|-|Tl|]. A standard theory of vortex dynamics has, however, not been 
established yet. For example, the existence of the Magnus force in superconductors is still 
an open question. In fact, Galilean invariance, which is used to derive the Magnus force in 
neutral superfluids, is unclear in charged superfluids because of the presence of a metallic 
lattice §. 

It is worth mentioning here that classical hydrodynamics provides an explanation 
of two main properties of vortices in a perfect fluid; (i) Helmholtz's theorem and Kelvin's 
theorem lead to a vortex motion such that a vortex-line moves with the background fluid, 
and (ii) Bernoulli's theorem leads to the Magnus force acting on a vortex-line at rest. The 
nonlinear Schrodinger approach || , which is one of the successful quantum descriptions for 
superconductors and superfluids, contains similar hydrodynamic equations: When the wave 
function is decomposed into the phase and the amplitude, ip = y/pe 9 , the field equation for 
9 becomes the continuity equation and the field equation for p gives Bernoulli's theorem. 
Since this theory contains Bernoulli's theorem, it should lead to the Magnus force. In 
fact, there exists an attempt to calculate the expectation value of the Hall (transverse) 
conductivity and the longitudinal conductivity due to vortex motions by solving the time- 
dependent nonlinear Schrodinger equation numerically |T(| . However, it is a non-trivial task 
to extract the qualitative features of vortex dynamics, in particular the Magnus force, from 
such simulations. 

In this paper, we develop an alternative approach to vortex dynamics, "the topological 
Landau- Ginzburg theory" which was proposed in a previous paper by two of us with Ao 
and Thouless ||. This is a field theoretical approach whose field equations include the 
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conventional nonlinear Schrodinger equation. In this approach, by using Helmholtz's and 
Kelvin's theorems as guiding principles, topological terms and collective coordinates for the 
vortex center are introduced. The hydrodynamical vortex law, "X = v(X)" , is naturally 
obtained by taking a variation of our "topological action" . Furthermore, the vortex motion 
can be extracted even after taking into account the phonon and photon effects. 

We will show that the topological term, u e flupX b flu fpx", and the source term in our action 
give rise to the Magnus force. This is consistent with the recent work by Ao and Thouless in 
which it is shown that the Berry phase and the dynamical phase in the quantum mechanical 
description of vortex system give rise to the Magnus force || . In our approach, the Magnus 
force appears naturally in both superconductors and superfluids, which is consistent with 
the expectation before M and with the recent results JlQUlI|] in different approaches. 



Our topological action also provides us with vortex-phonon couplings, where the phonon 
here is the density fluctuation of the order parameter at zero temperature, namely the 
zero sound wave. In actual superfluids and superconductors, there are situations where the 
vortices are trapped by pinning potentials instead of moving with the background flow. If 
the background flow becomes strong enough, depinning occurs and vortices start to move. 
A candidate of the driving force to overcome the pinning at zero temperature is the Magnus 
force and the dynamics of depinning is dictated by the inertial mass of a vortex. An 



alternative candidate is quantum tunneling | 17fl , but we do not consider it in this paper 



We will calculate the inertial mass as well as the vortex energy in neutral and charged 
superfluids by taking into account the phonon and photon contributions. Also evaluated 
are the core contributions to the inertial mass and the vortex energy in charged superfluids 
where the size of a vortex-core is generally larger than the atomic scale. The obtained 
inertial mass is applicable to a wider range of the parameter region than that obtained 
in the phenomenological time-dependent Landau- Ginzburg theory supplemented with the 
Fermi liquid theory ||T8 |. 



The organization of this paper is as follows. In section II, a close analogy between the 



Lorentz force in electrodynamics and the Magnus force in hydrodynamics is discussed using 
our topological action for the superfluid without vortex-phonon interaction. In section III, 
the effective vortex dynamics in a superfluid are studied by taking into account the phonon 
interaction. In particular, the inertial mass and vortex energy are calculated. In section 
IV, vortex dynamics in a superconductor are examined by taking into account the photon 
interaction and the density fluctuations. The inertial mass and vortex energy in conventional 
and high T c superconductors are then evaluated. Section V is devoted to summary and 
concluding remarks. 



II. THE MAGNUS FORCE 



In this section, by using the topological Landau-Ginzburg theory without phonon/photon 
fluctuations, we will demonstrate that the origin of the Magnus force in superfluids is quite 
analogous to that of the Lorentz force in electrodynamics. 

The action of topological Landau-Ginzburg theory for the neutral superfluid reads 
S = d 4 x 



1 2 

ip* (ihdo + ha ) if) - —\(ihdi + ha^ipl 2 - g (\^\ 2 - Po) 



2m 
2m 



J^(x) = J2 la J drda a d° <* (4) (* ~ r)) , (2) 

a=l 

where 7 a = 27rn a fr/m with integer n a and 27rh/m denotes vorticity unit. Notations follow 
from our previous paper 0. a M is a vector potential representing vortex singularities, b^ v 
is a rank-two anti-symmetric tensor potential, J^ u is the vorticity current, and U P i n (X) is a 
pinning potential. We examine vortex dynamics by taking a variation of S with respect to 
the vortex coordinate X. 

There is a one-to-one correspondence of our theory to electrodynamics where a point-like 
charged particle has a current 



U pin {X), (1) 



r d X^ 

r{x) = e / dr^— 5 (4) (x - X(r)). (3) 
J dr 

Variation of the source term — / d^xA^J^ in electrodynamics with respect to the particle 
coordinate X(t) gives the Lorentz force 

5 r 4 , v .dX"(t) 



SX^t) 



J d'xA.r = -eF t ,(X)^p- (4) 
e(£(X) + X x Bpf))* 



where the proper-time variable r is chosen to be r = t = X°. 

Analogously a variation of the source term / d A xbJ in (|l|) with respect to the vortex 
coordinate gives 

5 r w. , „ , ax"!* 4 ax""] 



/nX u* fix "1 
JxbrJ"' = -yH^X) — — , (5) 



where -ff^p = 9^6^ + + d p b^ v and r = t = X°. On the other hand, the field equation 
for a M obtained from (|]) gives a relation between the 6-field and the hydrodynamical current, 
J^{x) = (mp(x),mp(x)v' l (x)); 

Hfiup "^-{ivpXiJ • (6) 

The 6-field is the antisymmetric rank-two tensor potential whose exterior derivative is the 
observable hydrodynamical current. This property is analogous to the vector potential 
whose exterior derivative is electromagnetic field. Note that H ifJLV (X) in eq.(||) should be 
considered as an external field acting on the vortex at point X, which is analogous to the 
electromagnetic field acting on the charged particle in eq.(|J). Inserting eq.@ into eq.([5D 
leads to the conventional Magnus force 

x f flX^ ( 3X V \ 

Sm/^ = - ww-s- (7-^) 

= -mp(X)[(X - v(X)) x oj(X)} 1 

= ^MagnusiX), (7) 



where we have used the fact that (7^-) is the vorticity vector, uj = V x v. The Magnus 
force, which is analogous to the Lorentz force in electrodynamics, is thus the fundamental 
force in the hydrodynamic theory. 

If there are no pinning potentials, the variation of S at the "classical" level 
dS/dX l (t,a) = gives an equation of a vortex F l Magnus (X) = — mp(X)[(X — v(X)) x 
u)(X)] i = 0. Analogous to the case of the electrodynamics, J^{X) and H^ up (X) take values 
at point X without the contribution of a vortex sitting at X in the weak current limit; 
namely p(X) 7^ and v(X) 7^ 0. Thus one arrives at the conventional classical hydrody- 
namical relation X = v(X) which means that the vortex moves with the same velocity 
as the background flow. Even if one takes into account the full vortex contributions to 
J^{X), as far as one considers the finite size of a vortex core, it is possible to derive similar 
hydrodynamical relation in an averaged sense: 

In fact, if one uses a distribution with finite core size instead of the ^-function distribution 
given in (|^), the vortex centers X a and the current J~^{X) must be defined by averaging 
over the distribution as denoted in by (•) in @. For "classical hydrodynamical" cases (|§) 
further reduces to (X) = (v(X)) since (pv(X)) = (p(X))(v(X)) holds. One can imagine 
easily this situation for a case such as an isolated straight vortex line along the ^-direction 
and a constant background flows in the x-y plane. 

The above discussions are valid not only for neutral superfluids but also for charged 
superfluids since the bJ coupling is universal in both cases. Namely the Magnus force is a 
fundamental force in both neutral and charged superfluids. 

Up to this point, we have not taken into account the interactions of phonons or photons 
with vortices. They actually induce a kinietic term |m e //X 2 with m e // being the inertial 
mass. This means that, beyond the "classical" level of S, the acceleration term is induced 
from the equation of motion of X, i.e., m e ffX 2 = FMagnus(X) + F pin (X). This will be 



discussed in detail in later sections. 



III. EFFECTIVE ACTION FOR VORTICES IN NEUTRAL SUPERFLUID 



Now let us calculate the effective action for vortices in a neutral superffuid by taking into 
account the phonon interaction at zero temperature. We begin with the topological Landau- 



Ginzburg theory (|T|) and @ and define the phase and amplitude variables as ip = J p{x)e ielyX \ 
Thus the action 5* reads 



S= <Px 



hp (-d 9 + oo) - p- (W + a) 2 - ^-(Vp) 2 - g(p - po) 2 
2m timp 



h 

2m pa 



(9) 



If there is a uniform background flow Vb g =constant and vortex excitations do not exist, the 
stationary state of the superfluid is characterized as 







, P = Po 



9 = 9, 



m 



rn 



(10) 



In order to take into account the phonon fluctuation, we expand the amplitude and phase 
variables around the above stationary solution, p(x) = po + Sp, 9(x) = 9b g + 59. Then the 
full functional integral and the effective action S e ff(X) can be written as 



Z = J V[5p,59,a fl ,b f , u }e iS{p ' d ' a ^ b ^' x) = e iS ^f {x \ 



(11) 



where S(p, 9, a M , b^ v \ X) — J d 4 xC — U pin (X) with 

(ha 



£ = £ st + 5p 



— hvb g ■ a — - — a 



2m 



-{5p 59)G- 1 



rn 



-5pa-V59 + 0(5 3 ) + (bf, bJ) 



Cst = Po ( fta - hv bg ■ a - —a 



(12) 
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Here C s t is a stationary part of the action and G 1 is a Hermitian matrix denned as 



G- 1 



-2s + i^ A ) (-fido-fivtg-V) 



\ 



(hd + hv bg -V) ^-A ) 



(13) 



where we used the Coulomb gauge for the a-field and divergence free property of the back- 
ground flow. G is nothing but the Green's function for the fluctuations and has poles at 



2 \ 2 



to = v bg ■ k ± e(fc), with e (k) = c s k + 



2m) 



(14) 



with e(k) being the Bogoliubov spectrum and c s being the zero sound velocity 



m 



(15) 



For small momentum k, this reduces to the massless phonon mode propagating with sound 
velocity under the background flow. After the Gaussian integration with respect to 5p and 
SO, the effective lagrangian becomes 



2m 



•A 



a - v bg ■ a 



2m 



-a~ 



Xd + v bg -V) 2 + e< 



a Q - v bq ■ a 



2m 



-a 



(16) 



For interactions between vortices, only the phonon modes with small momentum are im- 
portant under normal circumstances, since the vortex motion and the background flow are 
much slower than the zero sound velocity; X « c s , v bg « c s . Therefore, from the second 
equality to the last one in (|T6|), (d + v bg ■ V) 2 is neglected and e 2 ~ — c 2 A is taken (the 
adiabatic approximation). 
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The integration of the b^ iv field is straightforward and the following constraint on is 



obtained: 



m 



f P x — d p a\ — d\a p — —Spx^J^ 



(17) 



This constraint (|TT| ) can be solved in the Coulomb gauge as 

a = ^(m/2h)e ijk d i J jk 



(18) 



A 



(m/h)e ijk d j J' 



OA- 



In the following, we will mostly focus our attention on the vortices in two spatial dimensions 
where vortices lie along the z-direction to make the argument as simple as possible. In this 



case, the a M field is written as 



a = - Ea n a e z ■ [X a x V In \x - X a \ 
a = J2a n ae z x Vln \x - X a \ 



(19) 



We will now insert the solutions (|T9|) into (|16[) to get an effective action S e ff(X) written 
in the vortex coordinate alone. When doing this, it is convenient to classify terms by the 
number of a^-fields. The a M field contains the vortex singularities, so the number of a M fields 
represents the number of the interacting vortices; for example, a bilinear term of a M field 
contains an interaction of two vortices. Let us first evaluate the term linear in a M field. 
Inserting the expressions ( |T9"D into the linear term in ([!(]), one gets 

J d 2 xp h(a -Vbg-a) = — J d 2 xp hJ2n a e z ■ (X a - v bg ) x Vln \x - X a \ 



mp 



lae z ■ (X a - v bg ) x (X a - Vbgt) 



(20) 



where one should be careful to perform spatial integration under the background flow. In 
fact, there is a constant ambiguity in the integral, / d 2 xVln\x — X\ = —tt(X + c). c 
can be fixed as —v bg t by requiring a correct boundary condition, i.e., the boundary of the 



S 



superfluid should also move with the background flow velocity. From (pO]), one sees that the 
combination a — v bg ■ a has the Galilean invariance. 



The obtained expression ( p0[ ) is a generalization of the Berry phase term [20| to the 
system with background flow. It is easy to show that taking a variation of this term with 
respect to the vortex coordinates leads to the Magnus force. 

Secondly, let us evaluate bilinear terms in the a M field. In order to evaluate them we will 
use the following formula: 

J d 2 xdi In \x — X a \dj In \x — X b \ 

-n{^H\X ab \/R) + {Xl b Xi b /X 2 ab ~^/2)} , ioiX a ^X b 
-n5 lj \n(d/R) , forX a = X b 



(21) 



where d is an ultraviolet cutoff which is an atomic scale, R is an infrared cutoff and X ab = 
X a — X b . Appendix A contains a derivation of the above formula by regulating In |x| by 
In \/x 2 + e 2 with d = ee 1 ^ 2 . It is natural to choose R to be the container size. In a charged 
superfluid, the penetration depth plays the role of the infrared cutoff R, since the charge 
screening makes the interaction region finite. 

Using the formula fl2"ID , / a 2 term in (|TJ) and fll6D becomes 



2m 



f d 2 xa 2 = (X>X)K n & / d 2 xV\n\x- X a \-V\n\x- X b \ 

•J 2iTTX i /i <J 



a=b a=£b 

-2 



p h 7T rV , 2 . JJ , r R . \X ab \ 



m " { ?^ ln 7 + 5^ (ln 7- ln ^ )} 

x ^ mpo \X ab \ 
-Eq + 2^ ^7a7b In —j- ■ (22) 



Here Eq is the static energy for a system with many vortices, 

£» = ^(I» 2 lnf . (23) 

The obtained expression Eq for an isolated vortex coincides with Feynman's result [2T|. If 



the container size is infinite, Eq for a single vortex is divergent, while Eq vanishes for a 



system with zero total vorticity. For example, a vortex — anti-vortex pair has zero static 
energy and has only the logarithmic interaction in (p^). 

Among the self-interaction terms, there is a term which is interpreted as the kinetic term; 
p h 2 



/^ X |^f (ao "^' a)2 

p h 2 



K J2 + J2)n a n b (X a - v bg ) l (X b - v bg ) j e ik e jl J d 2 xd k In \x - X a \d t In \x - X b \ 

s a=b a^-b 
■ l m effCE—( X *- V b 9 )) 2 

-?~\^llalb{{X a - v bg ) ■ (X b - v bg )(\n^j^- + ~) - (X a - v bg ) ■ u ab (X b - v bg ) ■ u ab ] 

°^ C s a+b d 1 

(24) 



where u ab = X ab /\X ab \. Here we define 

X cv = J2-(X a -v bg t) , 7 o = — (25) 
a 7o m 

and call this "the center of vorticity" of the system. This is analogous to the center of mass 
except that vorticity takes both positive and negative value. Then the first term of eq.(p4"D 
is interpreted as a kinetic term for Xqv with the following inertial mass 

m ^ = 4^ 7oln 7 • (26) 

If phonon fluctuations are neglected (i.e., in the limit c s — *> 0), "the center of vorticity" 
J2laX a / , ~fQ is conserved since it is proportional to the conserved total momentum P = 
J2p a = —{ m Po/2) e z x J2laX a . If there exist phonons, P has an extra term from the 
phonon momentum and J2laX a is not conserved anymore. Instead, it gets mobility with 
an inertial mass m e ff given in (pG|). 

The magnitude of the vortex energy Eq in (p3|) and the inertial mass m e ff in (p6|), are 



consistent with those in ref. |18j and ref. |T9[ obtained by different approaches, although 



the vorticity dependence was not considered in these references. Ref. [jn| is based on a 
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phenomenological time-dependent Landau- Ginzburg theory supplemented with the Fermi 
liquid theory, while ref. [19| is based on the functional integral approach at finite temperature. 
It is, however, unclear whether our result is consistent with ref. || which discusses the 
behavior of m e // for large R. 

Lastly we evaluate the vortex interactions. The leading term of two-vortex interaction 
is given by the second term of the right hand side of (|22|). Next to leading terms such as 
the second term in ([24|) are suppressed by the factor v\ g j(? s which is of order 10~ 4 ~ 1CT 8 
for superfluids He 4 , 10 -2 for conventional superconductors and 10~ 6 for high T c supercon- 
ductors. The suppression factor X /c 2 is even smaller than v 2 g / c 2 because we consider the 
adiabatic motion of vortices. Three and four vortex interactions are generated from third or- 
der and fourth order terms in a u . However, they can be neglected in a dilute vortex system. 
Furthermore these terms are suppressed by 1/c 2 and H. Third and fourth order terms also 
contain contributions to the vortex self-energy and the two-vortex interactions, but they are 
suppressed by h compared to leading terms in ( p2|) which are already of 0(h 2 ). 

The resultant effective action up to 0(h 2 ) and 0((l/c°)) is summarized as 

Seff(X) = J dt[L Berry + L ve + L 

kin -^int ] — U p i n (X), 

(27) 

L B erry = ^ Ea la^z ' (X a - V bg ) X (X a - V bg t) 

L ve = -E , with E = (mp /4ir)(E a la) 2 MR/d) 

< 

Lkin = \m efS X cv , with m eff = {mp Q / Anc 2 s )^lhi{Rl d) 

Lint = T,a^b( m Po/^halb^(\X ab \/d), 

where L Berry, L ve , Lkin and Li nt are the Berry phase term, the vortex energy, the vortex 
kinetic term and the interaction term, respectively. 
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The equation of motion for an isolated vortex with n a — 1 at zero temperature is deter- 
mined by varying the effective action which is the sum of the acceleration term, the Magnus 
force, and the pinning force; 

m e ffX = -mpo (X - v b Jj x u> + F pin 

F Magnus ~\~ F p ; m . (28) 

This equation gives the following scenario for vortex pinning and depinning. Consider a 
vortex pinned by the pinning potential U pin (X). The vortex feels the Magnus force as 
long as it does not move with the background flow v bg . If the Magnus force is strong 
enough to overcome the pinning potential, the vortex starts to move (depinning) with the 
accerelation dictated by the inertial mass m e //. As the vortex velocity gets closer to the 
background velocity v bg , the Magnus force gets smaller and reduces eventually to zero. The 
final stationary situation is then described by the hydrodynamical law for a perfect fluid; 
X = v(X). The phonon interaction is irrelevant at the final stage, but it is essential at the 
depinning stage. 



IV. EFFECTIVE ACTION FOR VORTICES IN CHARGED SUPERFLUID 



Now let us consider vortex dynamics in type II superconductors. We propose the follow- 
ing topological Landau- Ginzburg theory for vortices in superconductors; 



J d A x i)* ( 



i | ihd + ha - ^.1„) V - 2~ 



ihdi + hdi — -Ai) r 



-9 



Upi n (yX^ , 



-Po)' 



(29) 



J^x) = J2jajdTda ° d ° ° d ° 8W ( X - Xa (a,r)) 



(30) 



where q = 2e and m = 2m e with e and m e being the electron's charge and mass, A^ being 
the electromagnetic potential and F^ v = d il A v — d u A^. p iat is the electric charge of the 



12 



background lattice. In the stationary situation without vortices, the total electric charge 
vanishes locally; p iat + (q/c)p = 0. 

The Meissner effect and the quantization of magnetic flux can be easily checked by the 
variation of the action (|29|). The variation with respect to A 1 gives 



d^F^ + J l = , J e = -£l(-hVe -ha+ -A) . (31) 

mc c 



Taking curl of (|31~1) leads to 



:d + f )b = ^- < 32) 



where the penetration depth A is defined as 

A = ^mc 2 /p q 2 . (33) 

Eq. (^) represents the Meissner effect that the magnetic field has a penetration depth. 
Next let us integrate eq. (|32D in the region D surrounded by a closed path. 

ch f , ch . 
— / uj ■ an = — 
q Jd q 



f cti f cfl 

/ (A 2 D + l)B -dn = — u-dn= —2nn = (p n . (34) 
Jd q Jd q 



If we choose the integration region D larger than the penetration depth around the vortex, 
OB can be neglected because there are neither a magnetic field nor an electric current. 
Therefore the whole magnetic flux coincides with the integer multiple of unit magnetic flux 
4>o = (ch/q)2ir. 

The Hall (or transverse) voltage and the longitudinal voltage are produced depending on 
the direction of the vortex motion, which follow from the Maxwell's equations: If a magnetic 
flux moves, the magnetic field changes in time only through the flux motion and an electric 
field is produced by 

V x E(x) = ~B(x) = ~ c {X ■ V)B{x) = ~V x (X x B(x)) 

^E = --X x B . (35) 
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Now let us examine vortex dynamics in a superconductor by taking the variation of our 
effective action ( p^|) of the "classical" level. Since X dependence arises only through the 
source term as with the superfluid vortex system, the obtained force acting on vortices in a 
superconductor is also the Magnus force (]?]). The Magnus force ([?[) is rewritten in terms of 
the magnetic field through eq.( 



- mp (X - v) x u) = -p (v - X) x (1 + X 2 B)B , (36) 

c 

where the local fluid velocity is related to the electric current by v = (c/qpo)J e . When the 
vortex is almost at rest (X ~ 0) and the photon momentum is smaller than 1/A (X 2 OB ~ 0), 
the right hand side of (|36|) reduces to the Lorentz force, (l/c)J e x B |fj]. 



Now let us evaluate the effective action for vortices in charged superfluids by taking 
into account the phonon and photon fluctuations. Analogous to the neutral superfluid, a 
stationary solution without vortex singularities reads 



= , p = po , — VQ bg A bg = v bg and 6 bg = --prrvL 

m mc 2n y 



d^F^bg = ~Po + Plat = , d"F^ bg = -^povi g . (37) 
c c 

Deep inside bulk superconductors in three dimensions, neither the magnetic field nor the 
electric current present. On the other hand, in the surface region of a bulk superconductor 
and in a thin superconducting film, the electric current is non- vanishing due to the effect of 
the boundary conditions. These features are a consequence of the above equations. In the 
following, we will mostly consider a thin film placed in the x — y plane and assume that the 
transport current (the background velocity) is constant in the x — y direction. 



By expanding p(x), 9(x) and A^x) in ( |2D| ) around the above solution (^), we get an 
effective action 

Z = J V[6p, 59, 5A^, a M , ^^(aM^.V;*) = e <s.//W ; ( 38 ) 

where 
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C = C st + 5p 



1 



nao — nv bg ■ a a H 5 A ■ — a 



+^(5p 59 5A 5A)G- 1 



5p 
59 
5A 

\ 6A J 



m 



5pha ■ — (V59 - ^-5 A) + 0(5 3 ) + (bf, bJ) , 
m cri 



£ s t = Po I ha - hv bg ■ a - —a I - -F bg ^ u F{; g , 



(39) 



with Hermitian matrix G 1 



G 



-i 



(~ 2 9 + 4^ A) -h(d + v bg ■ V) -q/c \v bg 



h(d + v b9 -V) 



-q/c 



V 







-A 



-(□ 



(40) 



Here the Coulomb gauge for the a-field and <5v4-field are used. For zero background velocity, 
G has poles at 



(41) 



where e is the Bogoliubov spectrum (|14]). The massless Goldstone mode has been absorbed 
into massive mode. In order to take into account density fluctuations and the photon effect, 
we integrate out these fluctuations; 



C eff (a^,b^;X) « £ 



St 



"2^ ( %0 ~ Vb9 ' a) - L a J (d + v bg .Vy + e* + l/\i 



n 2 



h(a - v bg ■ a) - —a 



2m 



Po^ 2 1/A 2 



-a 



2m D + l/A : 



r a + 0(«7c2) 



(42) 
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Suppose that a thin superconducting sample is set in the x — y plane and the magnetic 
field is applied in the z-direction. Vortices lie along the z direction. Suppose also that the 
transport current J e flows in the x — y plane. By carrying out the integration of the 6-field 
and using the solutions (|19l), the first term of the effective action ( f42|) becomes 

p h 2 



2m 



« (tr) ^ + T,>an b (x a - v bg )\x b - W 

\ s I a= b a=£h 

Jd 2 xd k ln\ X -X a \ dl ln\ X -X b \ 



A x - l/(Ac s 

= E ~ m efJ,a( X a ~ V bg f 
a z 

-(*--^) ■ ( a- i/^).^" ) • (43) 

Here 27r-K ( a; ) = i^ 2 H (ix) with being a Hankel function of the first kind. The time 
derivative (c?q — v ■ V) 2 is neglected because X and i; are much smaller than c s , and the 
Bogoliubov spectrum e 2 is approximated by — c 2 A in the adiabatic approximation. In the 
last line of (pE3| ) the interaction part is neglected, because |X a b|c >> Ac,. 

One should note the qualitative difference between the kinetic term in the superfluid 
case ( |27|) and that in the superconductor case (|43|) . "Center of vorticity" enters into the 



kinietic term in the former, while each vortex has separate contributions to the kinetic term 
in the latter. In fact, in the superconductor, charge screening makes a vortex-size finite and 
each vortex can carry individual mobility and energy. The inertial mass of a single vortex 
is estimated using the previous regularization with the coherence length £ as an ultraviolet 
cut off, 

i fd>oV 

m eff,c 



1 ( 00 \ 
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So far, we have been adopting the vorticity distribution with the ^-function form (|30|), for 
simplicity. In this approximation, the ultraviolet cutoff should be the coherence length £. 
For neutral superfluids, £ is nearly equal to the atomic scale d which is the cutoff of the 
theory, so the 5-function vorticity is a good approximation. For superconductors, however, 
£ is larger than the atomic scale, so the vortex-core contribution could be non- negligible. In 
order to estimate the core contribution, let us separate the integral of the left hand side of 
(|jD as J r =£ d 2 x — > ( Jr = d + J* r =f ) d 2 x: The second term, which is a contribution from outside 
the core, is given by fl43|) and (|44]) . The first term from inside the core, can be evaluated 
in the following way. In the vorticity tensor, the 5-function is replaced by the regularized 
function, for example 

S M( x -X)^6(x -X )6(x 3 -X 3 )(-^j e -(^#) 2 , (45) 

which is introduced in |TjJ. So the vorticity vector u) has a finite peak at X with width £. 
By solving the constraint fll7D , one can show that the velocity around a vortex, a, has a 
maximum value at the edge of a core and decreases toward a vortex center, then vanishes at 
the vortex center. So we replace Vlnr in by (r/7r£ 2 )P(r/£) where P is a solution of (|I~7|); 
V • (rP(r/£)/7r£ 2 ) = e~^ x ~ / '(vr£ 2 ). One can also show that has linear dependence 
in r near the vortex center with P(0) = 1. Then the core contribution becomes 

2\Xc 8 ) a Jr=d ne " VA-l/(Ac a )V 7r£ 2 

-5te)'"Wfs:^ (46) 



where 



7T 2 J A 



a(A, n) = ± [' d 2 f rP(f) ■ (j^J rP(r) , (47) 
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with r = r/£. For a small core, £ ~ d, a is almost zero, since a(A ;/ u) »s |(1 — A)[rrP(r) ■ 
^ ^-P(^)]|f=(i-A)/2- For a large core, £ >> d, a becomes a constant of order 1, since 



A 



A-(£e/Ac s ) 



takes 1 for A >> (£c/Ac s ) , constant of order 1 for A m (£c/Ac s ) , and for 
A << (£c/Ac s ) 2 . The resultant effective mass including with the contribution from outside 
the core is 

<>*=U&)'* K *& + 4& ■ (48) 



Our result is consistent with that in ref. [Tj| for £c << Ac s . However, our formula is 



not limited to this parameter region, while that of is only valid in this region. For 
conventional superconductors, we have the Landau parameter A/£ ~ 10 and a sound velocity 
c s ~ 10 _4 c, while for high T c superconductors we have A/£ ~ 10 3 and c s ~ 10 _3 c. Thus the 
core of vortices in conventional superconductors is relatively large, d/£ ~ 10~ 4 , a ~ 10 -6 and 
K ~ 10~ 400 . On the other hand, that in high T c superconductors is small, c?/£ ~ 10° 1 , 
a ~ and i^o ~ 0.4. Therefore, the effective vortex mass in a conventional superconductor 
mainly comes from inside a core , while that of high T c superconductor mainly comes from 
outside the core. 

Now let us return to the evaluation of the a 2 term in S e ff(a^X). The second term of 
the effective action (f42|), when added to the term in the stationary action, £ st , becomes 



2m J \ □ + 1/X\ 

~ -TT 1^1 (Y^ + y^)n a n b - — ^ X — [ d 2 xdi In Ice — XJdjln \x — X b \ 

8n 2 \\) ^ ^ ;afc A x -l/A 2 7 41 a| 1 " 



In^^— , |X a6 | « A 

I -A ail 

^ e -|X ao |/A | Xa6 | A 

I -A- ab 



(49) 

where the vortex energy is estimated to be 
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— s(t)V»(|) 

The total vortex energy with the core contribution reads 

~rM n> I • ' <<x - (51) 

The first part of the resultant E^ 1 coincides with the free energy of a vortex line per unit 



length given by Abrikosov ||22|| . The core contribution is generally small; for the conventional 
superconductor case K ~ 3 and a ~ 0.1 and for the high T c superconductor case K ~ 1 
and a ~ 0. 

The vortex-vortex interaction in (pE9| ) has a natural form: At small distances, the same 
logarithmic force as the superfluid vortices acts between the vortices. For large distances, 
the force is exponentially suppressed because of the electromagnetic shielding in supercon- 
ductors. Other vortex interactions are suppressed by the factor v 2 /c 2 s . 

The resultant effective action for vortices in a charged superfluid up to 0(h 2 ) and 
0((l/c s )°) is given by 

S eff (X)= dt[ d 2 x C em + Y.{L Berry, a L ve a -\- -^fcjn.a} -^irat] Upi n , 
J J a 

(52) 
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LBerry,a = P^fin a e z ■ (X a - V bg ) X (X a - V bg t) 

< L ve , a =-E™ , with ^ = i L(0 o /A) 2 n2(K o (i) + a(|,I)) 
L km , a = \ml% a {X a -v bg y , with m *°; /ja = ^(0o/Ac s ) 2 ^(Xo(£) + a(f^)) 

where C em is the electromagnetic Lagrangian. 

The equation of motion for an isolated vortex is determined by varying the effective 
action which is the sum of the acceleration term, the Magnus force and the pinning force. 
A scenario for vortex motion is as follows: At first, a vortex in pinned, so that X is almost 
zero. If the transport current is large enough, the Magnus force becomes strong so that 
the vortex starts to move by overcoming the pinning potential. A longitudinal voltage is 
thus produced. The Magnus force acts in such a way that the vortex velocity gets closer to 
the transport velocity. Eventually the vortex moves along with the transport current. As a 
result, Hall (transverse) voltage is produced. 

V. SUMMARY AND CONCLUDING REMARKS 

We have presented a unified treatment of vortex dynamics at zero temperature by us- 
ing the topological Landau-Ginzburg theory. It is shown that the rank-two antisymmetric 
tensor potential, 6^^-field, is a fundamental field coupled to vortex lines, such that the local 
coupling "6^J MI/ " and the topological term "e M ^ pA 6 M ^/ ' p \ give rise to the Magnus force in 
both neutral and charged superfluids. Analogies with the electromagnetic theory have been 
shown; the rank-two antisymmetric tensor potential, the hydrodynamical current and the 
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Magnus force correspond to the vector potential, the electromagnetic field and the Lorentz 
force, respectively. 

Although we treat b^ u as a Lagrange multiplier in this paper, the action written in terms 
of b^y after integrating out the a M field in our theory ([I]) will have the same structure as 
the Kalb-Ramond type theory pl |. In this case, the 6^-field becomes dynamical. This 



may clarify the relation between our theory and that in (T^j23| where b^ is treated as a 

P 

fiuX 



dynamical field with a H^ uX term. 



The phonon effect in a neutral superfluid was calculated and our main result is the 
effective action given by fl27l) . In our derivation, u p Q h(a Q — Vb g ■ a)" becomes the Berry 
phase term after integrating out the 6^-field and it leads naturally to the Magnus force. 
In deriving the Berry phase term and the Magnus force, a first order formalism with single 
time-derivative is essential: Phenomenological Landau- Ginzburg theory having second order 



time-derivative (such as that in [18|) cannot lead to the Berry phase term. 



We provide a systematic method to regulate vortex self-interactions such as the vortex 
energy and the inertial mass. It turned out that the vortex energy and the inertial mass 
in neutral superfluids can only be defined for the whole system and not for the individual 
vortices. In fact, the vortex energy is proportional to the square of the total vorticity. Also, 
the real dynamical degree of freedom is "the center of vorticity" of the system as is seen 
from the effective kinetic term. 

The obtained energy Eq for an isolated vortex coincides with Feynman's chemical po- 
tential. The obtained inertial mass is about mpo7o/ c s ( m i s the atomic mass, po is the 
condensation density, 70 is the unit vorticity and c s is the zero sound velocity). They satisfy 
the relation; m e ff = E /cj. which is consistent with the result in [HJ where phenomenolog- 



ical Landau- Ginzburg theory is used supplemented with the Fermi-liquid theory. We have 
clarified the origin of the inertial mass by dividing fields into the stational configuration and 
the zero sound wave: a vortex follows the classical hydrodynamical law, U X = v" , if the 
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zero sound wave is neglected. Once the zero sound wave is taken into account, a vortex gets 
an inertial mass and deviation from the classical hydrodynamics arises. Therefore the zero 
sound wave should play an essential role for the pinning and depinning phenomena and also 
for the quantum tunneling of vortices. 

An effective action with photon and density fluctuations in charged superfluids is given 
by (|52D . The fundamental force acting on vortices in superconductors is also the Magnus 
force. For a slowly moving magnetic flux (vortex), the magnus force reduces to " J e x B 11 
which is the origin of the Lorentz force between the magnetic flux and the transport current. 
Vortices get inertial mass through 5A and 5p fluctuations. In contrast to the superfluid 
case, the vortex energy and the inertial mass can be defined for individual vortices because 
of the charge screening in superconductors. Core contributions to the vortex energy and the 
inertial mass are also evaluated. The obtained inertial mass is about lO _6 (l/47r)(0o/^ c s) 2 
for conventional superconductor and O.4(l/47r)(0 o /Ac s ) 2 f° r high T c superconductor, where 
0o is the unit magnetic flux and A is the penetration depth. Thus the main contribution 
to the effective inertial mass for conventional superconductors comes from inside the core, 
while that for high T c superconductor comes from outside the core. In reference opposite 
results are obtained for high T c superconductors. The reason for this discrepancy comes from 
the unjustified extrapolation of the formula in ]18|] from the region £c << Xc s to the region 
£c ~ Xc s . We do not have such a problem, since our formula is valid for all parameter space. 

Finally, we mention that the inertial mass of a vortex may be determined experimentally 
by the real-time observation of the vortex motion using electron microscopy p5f . 
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APPENDIX: A 



In this appendix we derive the integral formula (pl|): 

d 2 xdi In | a; — X a \dj In \x — Xf,\ 

-7r[S ij {H\Xab\/ee 1 / 2 ) - HR/ee 1 ' 2 )} + (X* ab Xi b /\X ab \ 2 - 5^/2)} , forX a ^ X b 



TiS ij \n(R/ee 



1/2 ^ 



, forX a = X t 



(Al) 



In the left hand side of ( |Alj) there is an ultraviolet divergence at X a = Xj, so we regulate 
In |x| as In \/x 2 + e 2 with a small parameter e. For an infrared divergence, we introduce 
a large cutoff parameter R so that we neglect 1/R 2 terms. Using Feynman's parameter 
formula 



1 

ab 



dt 



[at + 6(1 - t)Y 



(A2) 



the left hand side of (|A1|) becomes 
d 2 x- 



(\x-X a \ 2 + e 2 )(\x-X b \ 2 + e 2 ) 

1 dt I d 2 x' {xH + X:b{t ~ m f +XJabt) 



-7r[^'{-i + ln-|- 
2 



[x' 2 + t(l - t)Xi h + s 
In 



1 & 
+ 



\ 4 ' |X a6 | 2 



1/2- 


Jl/4 + e 2 /\X ab \ 2 


1/2 + 


y/l/4 + eVlX^I 2 



1 



\x ab \ 2 



l/4 + £ 2 /|X a6 | 2 



In 



1/2- Jl/4 + £ 2 /lX- 12 



- ab 



l/2 + Jl/A + e 2 /\X ab \ 2 



)] 



(A3) 



It is easy to recover the right hand side of flAj]) for both X a X b with \X ab \ >> e and 
X a = X b . Also, the above formula (A3) may be useful for a case where two vortices are 
close to each other \X ab \ ~ s. 
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